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Abstract 

In this paper we present a differential variational inequality formulation of dynamic net- 
work user equilibrium with elastic travel demand. We discuss its qualitative properties 
and provide algorithms for and examples of its solution. 
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1 Introductory remarks 

This paper presents a differential variational inequality representation an elastic de- 



mand exten s ion of the fixed-demand dynamic traf fic assignment model original! 
Friesz et al.l (119931') and discussed s u bsequ ently by iFriesz and Mookherjed (|20O 



presented in 



Friesz et al. 



(l201ll . 120121 ): IFriesz and Meimandl toij ). As such, it is concerned with a specific type of 
dynamic traffic assignment known as dynamic user equilibrium (DUE) for which travel cost, 
including delay as well as early and late arrival penalties, are equilibrated and demand is 
determined endogenously. 

2 Some history 

In this section we review some of the fe w prior efforts to model D UE with elastic travel 
demand. This review is based in part on IFriesz and Meimandl (|2012l ). Most of the studies 
of DUE reported in the DTA literature are about dynamic user equilibrium with constant 
travel demand for each origin-destination pair. It is, of course, not g enerall y true that travel 
deman d is fixed, even for short time horizons. Arnott et al. ( 19931 ) and Yang and Huang 
(Il997h directly consider elastic travel demand; their work possesses a limited relationship to 
the analysis presented in this paper, for their work is concerned with a simple bottleneck 
instead of a nontrivial network, which is our focus. 
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Yang and Meng (1 19981 ) extend a simple bottleneck model to a general queuing network 
with known elastic demand functions for each origin-destination (OD) pair. IWie et al.l ((2002) 
study a version of dynamic user equilibrium with elastic demand, usin g a complementarity 
formulation that requires path delays to be expressible in closed form. ISzeto and Lol (|2004l ) 
study dynamic user equilibrium with elastic travel demand when network loading is based 
on the cell transmission model (CTM); their formulation is based on discr ete time and is 
expressed as a finite-dimensional variational inequality (VI). Han et al. ( 201ll ) study dynamic 
user equilibrium with elastic travel demand for a network whose traffic flows are also described 
by CTM. 

Although IFriesz et al. I (mil) 

show that analysis and computation of dynamic user equilib- 
rium with constant travel demand is tremendously simplified by stating it as a a differential 
variational inequality (DVI), they do not discuss how elastic demand may be ac commodated 



within a DVI framework. We know of no prior published work that extends the IFriesz et al. 



(mil) formulation to an elastic demand setting. That is, our demonstration in this paper 
that dynamic user equilibrium with elastic travel demand may be formulated as a differential 
variational inequality in continuous time is original and has not been previously reported. 

The DVI formulation of elastic demand DUE is not straightforward. In particular, the 
DVI presented herein has both infinite-dimensional and finite-dimensional terms; moreover, 
for any given origin-destination pair, inverse travel demand corresponding to a dynamic user 
equilibrium depends on the terminal value of a state variable representing cumulative depar- 
tures. 

The DVI formulation achieved in this paper is significant because it allows the still emer- 
gening theory of differential variational inequalities to be employed for the analysis and com- 
putation of solutions of the elastic-demand DUE problem when simultaneous departure time 
and route choice are within the purview of users, all of which constitutes a foundation problem 
within the field of dynamic traffic assignment. A good review of recent insights into abstract 
differential variational in equality theory, in c ludin g computational methods for solving such 
problems, is provided by iPang and Stewart (|2008T ). Also, differential variational inequalities 
involving the k ind of explicit, agent-specific control variables employed herein are presented 
FrieszJ (|2O10h . 



in 



3 Notation and essential background 



We assume the same general setup as IFriesz and Meimand ( 20121 ). In particular, the time 
interval of analysis is a single commuting period or "day" expressed as 

[t ,t f ] c 

where tf > to, and both to and tj are fixed. Here, as in all DUE modeling, the single most 
crucial ingredient is the path delay operator, which provides the delay on any path p per unit 
of flow departing from the origin of that path; it is denoted by 

D p (t,h) for aliped 

where V is the set of all paths employed by travelers, t denotes departure time, and h is a 
vector of departure rates. From these, we construct effective unit path delay operators ^ p (t, h) 
by adding the so-called schedule delay E[t + D p (t, h) — Ta\; that is 



tf p (t, h) = D p (t, h) + E[t + D p (t, h) - T A ] for all p € V 
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where T4 is the desired arrival time and T4 < tf. The function E {•) assesses a penalty 
whenever 

t + D p (t,h)^T A (3.1) 

since t + D p (t, h) is the clock time at which departing traffic arrives at the destination of path 
p € V. We stipulate that each 

tf p (-,/i) : [to,*/] — for aliped 
is measurable and strictly positive. We employ the obvious notation 

(*„(-, /i) :peV)e 

to express the complete vector of effective delay operators. 

It is now well known that path delay operators may be obtained from an embedded delay 
model, data combined with response surface methodology, or data combined with inverse 
modeling. Unfortunately, regardless of how derived, realistic path delay operators do not 
possess the desirable property of monotonicity; they may also be non-differentiable. However, 
we shall have more to say about path delays in a subsequent section of this paper. 

Transportation demand is assumed to be expressed as the following invertible function 

Qa(tf) = FiM 

for each origin-destination pair £ W, where W is the set of all origin-destination pairs and 
v is a concatenation of origin-destination minimum travel costs associated with € W. 
That is, we have that 

Vij € 3?+ 
v = ( Vij : € W) €^ W I 



Note that to say is a minimum travel cost means it is the minimum cost for all departure 
time choices and all route choices pertinent to origin-destination pair G W. Further 

note that Qij (tf) is the unknown cummulative travel demand between € W that must 
ultimately arrive by time tf. 

We will also find it convenient to form the complete vector of travel demands by concate- 
nating the origin-specific travel demands to obtain 

Q = (Q i ,(t / ):( i ,j)eW)eKl >v l 
F : Rfl_>Hfl 

The inverse demand function for every € W is 

= Qij [Q (t f )}, 

and we naturally define 

& = (&ij :(»j')eW)e#l 

Additionally, we will define the set Vij to be the set of paths that connect origin-destination 
pair (i, j) G W. We denote the space of square integrable functions for the real interval [to,tf] 
by L 2 [to,*/]- We stipulate that 

he(Ll[t ,tf]f l 
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We write the flow conservation constraints as 

rtf 

hp (t) dt = Qij {tf) for all G W (3.2) 

pePn Jt ° 

where (|3.2I) is comprised of Lebesgue integrals. In (|3.2h . we consider demand to be the cum- 
mulative travel that takes place by the end of the analysis horizon tf. We define the set of 
feasible flows for each v G 5ft by 

Ao = \ h > : h P ^ dt = Q*i (*/) for a11 G W f ^ ( L + 

Given h let us also define the essential infimum of effective travel delays 
Vij = ess inf h) : p £ Vy] for all (i, j) G W 

We rely on the following definition of dynamic user equilibrium: 

Definition 3.1. Dynamic user equilibrium. A vector of departure rates (path flows) h* G Ao 
is a dynamic user equilibrium if 

h* p (t) > 0,p G Pij =► [t, h* (t)} = Vij 

We denote this equilibrium by DUE Ao, [to,tf]). 

The meaning of Definition 13.11 is clear: positive departure rates at a particular time along a 
particular path must coincide with least effective travel delay. Furthermore 

^ p (t, h*) > Vij,p G Pij ^h* p = 

is an implication of Definition 13.11 



4 Dynamics 

We assume there are unknown terminal state variables Qij(tf), for all G W, which are the 
realized DUE travel demands. Moreover, for each origin-destination pair G W, inverse 
travel demand is expressed as 

v lj = e iJ [Q(t f )] (4.3) 

Thus, the Qij(tf), for all (i,j) G W, will be determined endogenously to the differential 
variational inequality presented subsequently in Section Such an approach contrasts to 
the approach employed by Friesz et al. (2010) to study fixed-demand DUE by making each 
Qij(tf) an a priori fixed constant. Accordingly, we introduce the following dynamics: 

^ = J2 h P <5y(*o) = for all (i,j) G W (4.4) 
As a consequence, we employ the following alternative form of the feasible set: 
A = 



h>0:^=J2h p (t) Qij(t o )=0 for all GW\ C {L 2 + [t ,t f ]y " (4.5) 

' pePi-j 



Note that the feasible set A in (14. 5h is expressed as a set of path flows since knowledge of h 
completely determines the demands that satisfy the inital value problem (|4.4p . 
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5 The Differential Variational Inequality 



Experience with differential games in continuous time suggests that an elastic demand dynamic 
user equilibrium is equivalent to the following variational inequality under suitable regularity 
conditions: find h* € A such that 



E 



'to 



%{t, h*)(h p - h* p )dt - J2 % iQ* (*/)] [Qh (</) - Qh (*/) 

(«)eW 



> for all h E A 



(5.6) 



where Q*(tf) is the terminal equilibrium demand and 0* [Q* (tf)] is the inverse demand func- 
tion evaluated at Q*(tf). We refer to differential variational inequality (|5.6p as DVI{^ , 0, to,tf). 
We now proceed constructively to show that (15, 6ft is equivalent to elastic demand DUE for 
appropriate regularity conditions. 

In particular we note that the differential variational variational inequality of interest may 
be written as 



E E 1* %{t,h*)h p dt- %iQ*(-tf)]Qij(tf) 



* E E f f y P (t,h*)h* p dt- Y %{Q*(tf)\Qi 



(tj)ew 

for all /i € A. Inequality ()5.7[) means that the solution h* £ A satisfies the optimal control 
problem 



min Ji 



o 



£ e&[Q* (*/)]<&(*/) + E E [ tf %(t,h*)h 



(5.8) 



subject to 



(i,i)eW 



^ for all (i,j)e>V 



Qij-Cto) = for all (i,j)eW 
> 



(5.9) 

(5.10) 
(5.11) 



This optimal control problem may not be used for computation, since it involves knowledge 
of the DVI solution. However, it may be used to express necessary and sufficient conditions 
for the solution of DVI(^f, 0, to, tf). In particular, the Hamiltonian for DVI{^, 0, to, tf) is 

H = E E^' /i *) /i p+ E E h p 

(i,j)eWp&Vij (ij')€W P&Vij 



(i 



E I E + 



A; 



where the adjoint equations are 



eft 



<9F 



for all (i,j) e W,p G V&t € [t ,t f ] 



(5.12) 



5 



with transversality conditions 



for all € W,p G P -,t G [t ,*/] (5.13) 



MV) = - t,Jjew JT.;,./ = -e?,- [Q* (*/)] 



It is clear from ([5TT2]) and ([5TT3]) that 

Ay (i) = — @*j [Q* (tf)] , a constant 

We note that the Hamiltoni an is l i near i n h and does not depend explicitly on the state 
variables. By Theorem 3.7 of Friesz ( 20ld ). the Mangasarian sufficiency theorem assures the 



minimum principle and associated necessary conditions are also sufficient. 

Since h is a control vector and must obey the minimum principle in for each instant 
of time, we enforce 

mini/ s.t. — h < 

h 

for which the Kuhn- Tucker conditions are 

%(t,h*) - % [Q* (t f )\ = Pp >0 for all G W,p G V%j , t G [t ,t f ] (5.14) 

where the p p are dual variables satisfying the complementary slackness conditions 

p p h p = for all G W,p G Vij , t G [t ,t f ] (5.15) 

From (|5.14p and (|5.15p we have immediately 

h* > 0, P EV ij ^%(t,h*) = @* j [Q*(t f )} 

%(t,h*) > e?.[Q*(t / )], P €7'i i =^/i; = o 

which are recognized as conditions describing a dynamic user equilibrium. 

We have noted above that (|5.14p is equivalent to a dynamic user equilibrium; thus, showing 
that (|5.14p corresponds to a solution of DVI(^f, 0, to,tf) will complete the demonstration that 
DVI(^>, ®,to,tf) is equivalent to a dynamic user equilibrium. In particular note that 

p p (h p - h*) = pphp - pph* = pphp > 

so that 

jttpfofc*) - Qtj iQ* (*/)]| ~h* p )>0 for all G W,p G V tJ ,t G [t Q ,t f ] (5.16) 

where h, h* G A. From (|5.16p we obtain 



%(t,h*)(h p - h* p ) - % [Q* (t f )} (hp -h* p )>0 for all (i,j) G W,p G V ih t G [t ,t f ] (5.17) 
From (|5.17|) . we obtain 
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Therefore 



W [ tf * p (t,h*)(h p -h* p )dt- £ %[Q*{t f )]f f\^9ll-^jl)dt>0, (5.19) 
peP J ^ (ij)ew J Jt ° 

from which , 6, to, tf) is obtained immediately, since Qij(to) = Q*j(to) = 0. 

The preceding analysis has established the validity of the following theorem: 

Theorem 5.1. Elastic demand dynamic user equilibrium is equivalent to a differential varia- 
tional inequality. Assume *S> p (-,h) : [fojif] — ► ^++ * s measurable and strictly positive for all 
p € V and all h £ A. Also assume that the elastic travel demand function is invertible, with 
inverse Qij (Q) for all £ W. A vector of departure rates (path flows) h* € A is a dynamic 
user equilibrium with associated demand Q*(tf) if and only if h* solves DVI($?,Q,to,tf). 



6 Additional qualitative properties 



We are presently working on existence results, which will be ready by the time of submission 
of the full paper. Those results will be wholly original. Also existence will be discussed 
from the point of view of the path delay operators that are assumed or derived from network 
loading consideratio ns. We will conduct analysis on several network flow mod els including 



the V ickrey model (jVickrevl . 119691; lHan et al 



2012a 



20121 ). the Link D elay Model iFriesz et al.l (|l993l): Irian et al.l (|2012dl ). and the kinematic wave 



bj), th e IWR -Lax model (jFriesz et al 



model (|Han et all l2012eh . We have analyzed the continuity of the effective delay operators in 



Hilb er t space . In particu lar , such an analytical pro perty has been est ablished for the Vickrey 
model ( Han et al. . 2012c ) and the link delay model Han et al. ( 2012dl ). We have also acquired 
the capability of analytically investigating the kinematic wave model on networks with vehicle 
spillback. In particular, existence, uniquenes s as well as well-p osedness of a network-based 
kinematic wave model has been established in lHan et all (|2012eh . 

These qualitative results on network loading models, when com bined with the general 
framework of proving DUE existence proposed by Han et aD ( 2012c ). will easily lead to the 
existence results of DUE with different network performance models. 



7 Computation 



The computation of solutions to (|5.6p is also under development and will be ready by the 
time of submission of the full paper. The algorithms and calculations will be wholly original. 
We have successfully created analytical DUE models supported by an e fficient DUE solution 
algorithm that has been extensively tested, see Friesz et al. ( 20121 . 2011 ). We expect to finish 
and demonstrate a DUE solution associated with the kinematic wave model with spillback for 
large-scale network by the submission of the full paper. 



8 Conclusion 

We have shown dynamic network user equilibrium based on simultaneous departure time and 
route choice in the presence of elastic travel demand may be formulated as a differential 
variational inequality (DVI). The above conclusions will be expanded to reflect the qualitative 
and numerical analyses of Sections [6] and El 
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